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Abstract 

Tiny violations of the Lorentz symmetry of relativity and the associated discrete CPT symmetry 
could emerge in a consistent theory of quantum gravity such as string theory. Recent evidence 
for linear polarization in gamma-ray bursts improves existing sensitivities to Lorentz and CPT 
violation involving photons by factors ranging from ten to a million. 
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Observations of photon behavior provide crucial probes of fundamental physics. Famous 
examples include the classic Michelson-Morley, Kennedy-Thorndike, and Ives-Stilwell ex- 
periments [1], which support the foundational Lorentz invariance of relativity. In recent 
years, a wide range of astrophysical, solar-system, and laboratory tests of Lorentz symme- 
try and its associated discrete CPT symmetry have achieved impressive sensitivities using 
photons (see Ref. [2] for a compilation). One motivation for these efforts is the prospect 
that tiny violations of these invariances could emerge in a consistent theory of quantum 
gravity such as string theory [3]. In this paper, we use recent measurements of linear po- 
larization in light from gamma-ray bursts (GRB) [4-6] to improve existing sensitivities to 
a variety of Lorentz- and CPT-violating effects by factors ranging from ten to a million. 
The key to the exceptional GRB sensitivity to Lorentz and CPT violation lies primarily in 
the extreme propagation distances during which tiny effects can accumulate, along with the 
comparatively high photon energies involved. 

At accessible energies, violations of Lorentz invariance are governed by the Standard- 
Model Extension (SME) [7], a comprehensive effective field theory containing both General 
Relativity and the Standard Model that provides a general theoretical framework for obser- 
vational studies. This theory also describes CPT violation in the context of realistic field 
theory [8]. The SME action is a sum of coordinate-invariant terms, including ones formed 
from Lorentz-violating operators contracted with controlling coefficients, and the mass di- 
mension d of each operator fixes the dimensionality of the corresponding coefficient [9]. In 
the photon sector, all gauge-invariant operators describing the propagation of light have 
been classified and enumerated for arbitrary d [10]. 

The SME predicts that light propagates in the presence of Lorentz and CPT violation 
as the superposition of two normal modes that may differ in speed and polarization. The 
dispersion relations connecting the photon energy E and momentum p for the two modes 
can be written in the compact but implicit form [7, 10] 

E=(l-q°± V(0 2 + (O 2 + (<T 3 ) 2 )P, (1) 

where the dimensionless quantities q a = q a (E,9,<f)) depend both on E and on the photon 
direction of propagation, which for an astrophysical point source is fixed by the source 
codeclination 9 = (90° — declination) and by the right ascension 0. For light propagating in 
vacuo, the quantities q a are linear combinations of the basic SME coefficients c|^ m , fc|^ m , 
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k(B)jm-> anc ^ k(v)jm-> where j, m are angular quantum numbers. For each d, many different 
SME coefficients control the behavior of light. A given point source at a fixed sky location can 
therefore access only a limited number of coefficient combinations. Consequently, multiple 
sources at different sky locations are required to constrain fully the coefficient space. 

Two major features, dispersion and birefringence, can be exploited to search for Lorentz 
violation in radiation from sources at cosmological distances. Dispersion is a characteristic of 
all SME operators with d ^ 4. Arrival-time differences in high-energy photons from sources 
such as GRB can be used to constrain the energy dependence in the group velocity. The 
SME framework shows that dispersion for operators with odd d is necessarily accompanied 
by birefringence, implying different speeds for the two normal modes. As a result, a wave 
packet not only disperses but gradually splits into two. Only the CPT-even operators with 
even d characterized by the coefficients c^- m give dispersion without birefringence. GRB 
constraints on Lorentz-violating dispersion for even dimensions d have been obtained for 
d = 6 and 8 [10-13]. 

Birefringence studies of astrophysical sources such as GRB offer extreme sensitivity to 
Lorentz and CPT violation. The primary signature of birefringence is a change in photon 
polarization due to propagation. This is governed by the phase difference of the eigenmodes 
developed during propagation, which increases with energy. While dispersion can uniquely 
constrain the coefficients c^ jm , for d > 4 birefringence at GRB energies typically offers 
many orders of magnitude better sensitivity to the coefficients k^^ m , fc[m jm , ^[yw because 
a comparable dispersion test would require a time resolution comparable to the tiny inverse 
photon frequency. 

In one study, evidence for polarization at the level of II > 35% in GRB 930131 and II > 
50% in GRB 960924 was extracted from data obtained by the Burst and Transient Source 
Experiment (BATSE) [14]. These results have been used to constrain Lorentz violation 
for d = 4,5,6,7,8,9 [10, 15]. Polarization as high as n = 96^% was identified in GRB 
041219A using instruments aboard the International Gamma-Ray Astrophysics Laboratory 
(INTEGRAL) [4]. This result was used by Stecker to place bounds on d = 5 coefficients on 
the order of 10" 34 GeV" 1 [16]. Another INTEGRAL analysis found similarly high degrees 
of polarization in GRB 041219A, placing bounds on the single isotropic d = 5 coefficient 
[17]. A recent analysis using polarization data for GRB 100826A (> 6%), GRB 110301A 
(> 31%), and GRB 110721A (> 35%) from the Gamma-ray Burst Polarimeter (GAP) on 
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z 


energy 


(M) 


GRB 930131 


0.1 [15] 


31 - 98 keV [14] 


(98°, 182°) [19] 


GRB 960924 


0.1 [15] 


31 - 98 keV [14] 


(87°, 37°) [19] 


GRB 041219A 


0.02 [17] 


100- 1000 keV [4] 


(27°, 6°) [20] 


GRB 100826A 


0.71 [18] 


70 - 300 keV [5] 


(112°, 279°) [21] 


GRB 110301A 


0.21 [18] 


70 - 300 keV [6] 


(61°, 229°) [22] 


GRB 110721A 


0.45 [18] 


70 - 300 keV [6] 


(129°, 333°) [23] 



TABLE I: GRB for which strong evidence of linear polarization exists. The second column gives the 
estimated lower limit on the red shift. The third column is the energy range in which polarization 
is observed. The last column gives the GRB codeclination 9 and right ascension <p. The first two 
GRB were previously studied [10, 15], while the others are the subject of the present work. 




FIG. 1: Sky map showing the locations of the six GRB listed in Table I. The map is centered at 
(0,0) = (90°, 0°). 

the Interplanetary Kite-craft Accelerated by Radiation of the Sun (IKAROS) [5, 6] also 
bounded the single isotropic d = 5 coefficient [18]. The basic features of all six GRB are 
summarized in Table I, and their celestial locations are displayed in Fig. I. 

In this work, we use the polarization reported for the four latest GRB to place improved 
bounds on direction-dependent combinations of SME coefficients for d = 4,5, 6, 7, 8, 9. For 
the analysis, it is useful to expand the quantities <, a in energy E [10, 11], 

, a {E,6,<P) = Y,E d -\ Wa {e,<P), (a = 0,1, 2, 3), (2) 

d 

where (9 , 0) are direction-dependent combinations of SME coefficients. The direction 
dependence can be displayed explicitly by further expansion using conventional spherical 
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coeff. 


CPT 


d 


j 


biref. 


general 


c (d) 

(I)jm 

k (d) 
(E)jm 

(B)im 

k (d) 
(V)jm 


+ 
+ 
+ 


4,6,8,... 
4,6,8,... 
4,6,8,... 
3,5,7,... 


0,l,...,d-2 
2,3,...,d-2 
2,3,...,d-2 
0,l,...,d-2 


/ 
/ 
/ 


isotropic 


(d) 
C (I)00 

l,(d) 
K (V)00 


+ 


4,6,8,... 
3,5,7,... 






/ 



TABLE II: Index ranges and properties of coefficients for Lorentz and CPT violation, where — j < 
m < j as usual. The isotropic limit is shown in the last two lines. 

harmonics Y jm (6,(f)) for ^°(9,<f>) and <?^ 3 (#,0), 

jm 

T (d)3 (M) = (3) 

jm 

and their cousins ±?Yj m (9, 4>) of spin- weight two for the combinations = q^ 1 =F iqW 2 , 

0) = ±2Y jm (9, 0) (k[% jm ± ifc«.J . (4) 

jm 

The basic SME coefficients - m , k^ jm , and fc[^L m characterize CPT-preserving Lorentz 
violation, while k^ jm also controls CPT violation. Table II summarizes some properties of 
these coefficients. 

Table II also contains information about the isotropic coefficients for which j = m = 0, 
which represent a popular restriction of the general framework. In this limit, each value 
of d has exactly one SME coefficient, which is nonbirefringent in the CPT-even case and 
birefringent in the CPT-odd case. The theoretical motivation for this restriction is open 
to doubt because isotropy can hold only in a single inertial frame, which cannot be an 
Earth-based frame and requires fine tuning to match the standard Sun-centered frame. 
However, the isotropic limit does offer an order-of-magnitude measure of the reach for a 
given source, and it also simplifies many equations. For example, the general expression 
for the defect in the group velocity relevant for dispersion studies can be written 5v gr = 
(d - 3)E d - 4 ( - ± in the CPT-even case and 5v gr = ±(d - 3)E d - 4 \^ 3 \ in 



the CPT-odd case, but in the isotropic limit one obtains instead the simpler expressions 
Svgr = -(d - 3)E d ~ 4 c\f )Q0 /V4^ in the CPT-even case and Sv gr = ±{d - 3)E d - 4 k { $ )00 / 
in the CPT-odd case. Note that Table II reveals there are no j = isotropic coefficients 
k(E)jm or ^(B)jm' s0 ^h e quantities are necessarily direction dependent. 

During birefringent propagation of light, the difference in phase speed between the two 
normal modes generates a relative phase shift and a corresponding change in the net polar- 
ization. This change can be visualized as a rotation of the Stokes vector s — (s 1 ,s 2 ,s 3 ) = 
(Q,U,V) about a rotation axis <f = (V,<; 2 ,^ 3 ) by an angle $ equal to the total relative 
phase shift [24] . If s lies along <f, then the light is in one of the two normal modes and 
the polarization remains unchanged. Note that in the CPT-odd case the two normal modes 
are circularly polarized, <f = (0,0,^ 3 ), while in the CPT-even case they are linearly polar- 
ized, <f = (V,^ 2 ,0). Also, isotropic birefringence occurs only in the CPT-odd case where 
^3 _ E d ~ 4 k^ Q0 / y/4ar and for which d is odd, so changes in polarization for even d necessarily 
depend on the source position. 

For given d, the rotation angle <3> depends on the difference in phase speed, Av = 
2E d ~ 3 \q( d ) a \, where a = 3 in the CPT-odd case and a = + in the CPT-even case. It can 
be found by integrating the accumulated phase over the propagation time, $ = J EAv dt. 
Incorporating the redshift gives 

$ = 2E d ~ 3 L {d) \q {d)a (9,<j))\, a=(3,+), (5) 

where = J**(l + zY'^H' 1 dz is the effective baseline for dimension d in terms of the 
source redshift z and the Hubble expansion rate H z at redshift z. 

For the CPT-odd case with a = 3, the Stokes vector rotates about the s 3 axis, shifting the 
linear polarization angle ip by Sip = $/2. A detailed analysis searching for CPT violation 
could take advantage of the E d ~ 3 dependence in Eq. (5). For example, assuming two photons 
of energies E\ and E 2 initially have the same polarization, the difference ip2 — i>\ in their 
polarization angles after traveling the effective baseline is given by 

M-t_ Jd-3\r(d) = S 0) k (V)jm- ( 6 ) 



(E d ~ 3 - E d ~ 3 )U d ) 



jm 



Using observations of GRB 041219A and assuming only the single isotropic coefficient with 
d — 5 [25], a search of this type was performed in Ref. [17], yielding the bound |fc(y) 00 | < 
3 x 1CT 33 GeV" 1 . 
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An alternative approach offering conservative constraints is to seek a significant de- 
gree of linear polarization within a given energy band. Differential rotations within the 
band smear the polarization and hence decrease the effective degree of linear polarization, 



Il e fj = ^ {s 1 } 2 + (s 2 ) 2 . In the CPT-odd case, this smearing produces an upper bound in 



the measured polarization II < y (cos$) 2 + (sin$) 2 , where the equality holds for an initial 
100% linear polarization at constant angle. The polarization smearing will be nearly com- 
plete unless the change in <£> over the energy band is less than 2n. Under the assumption 
of only the single isotropic coefficient with d = 5, a similar idea has been applied to ob- 
servations of GRB 100826A, GRB 110301A, and GRB 110721A [18], producing the bound 

Kvoool < 6 x 10 " 34 GeV_1 - 

In the general case, for any given odd d and including all coefficients for CPT violation, 

we obtain the conservative limit 



< * (7) 

\Et 3 - Et 3 \LW' { ) 

where E\ and E 2 are the edges of the energy band. This expression can be used to obtain 
substantially improved sensitivities to CPT violation from GRB 041219A, GRB 100826A, 
GRB 110301A, and GRB 110721A. We consider CPT-violating operators of dimensions d = 
5, 7, and 9, for which there are 16, 36, and 64 independent vacuum coefficients, respectively 
[10]. Each source generates a new constraint on a different combination of coefficients, listed 
in Table III. For completeness, we also list bounds in the isotropic limit. These results 
represent improvements in sensitivity over existing limits on violations of CPT symmetry 
ranging from tenfold to about a millionfold. 

The CPT-even case is more complicated because the normal modes are linearly polarized, 
which implies linearly polarized light from a distant source produced near one of the two 
polarizations s « ±(V, <^ 2 , 0) could propagate essentially unchanged. Moreover, even if both 
modes are involved, the change in polarization involves more than a simple rotation of linear 
polarization. Light that is initially linearly polarized, s — (s^s^O), becomes elliptically 
polarized as s rotates out of the s 1 -s 2 plane and in some cases may even become circularly 
polarized, s— (0,0, s 3 ). 

For light not produced in a normal mode and initially linearly polarized at angle tpo, 
let \I/ = ipo — ipb be the difference between ip and the polarization angle ipb for the faster 
of the two normal modes. As the Stokes vector rotates about <f, it traces out a cone with 





GRB 041219A 
(27°, 6°) 


GRB 100826A 
(112°, 279°) 


GRB 110301 A 
(61°, 229°) 


GRB 110721A 

(129°, 333°) 


|E im >W(M) fcg>,.J 


< 2 x 10" 34 GeV" 1 

< 2 x 1(T 28 GeV" 3 

< 2 x IO" 22 GeV" 5 


< 7 x 1(T 35 GeV" 1 

< 4 x ltr 28 GeV" 3 

< 2 x IO" 21 GeV" 5 


< 2 x IO" 34 GeV" 1 

< 2 x 10- 27 GeV" 3 

< 2 x 10- 20 GeV" 5 


< 1 x IO" 34 GeV" 1 

< 8 x IO" 28 GeV- 3 

< 5 x IO" 21 GeV- 5 


lfc (s) 1 

l K (V)Ool 
1,(7) 1 

r(v)ool 
r(v)ool 


< 8 x IO" 34 GeV" 1 

< 8 x 10~ 28 GeV -3 

< 8 x IO" 22 GeV" 5 


< 2 x 1CT 34 GeV -1 

< 1 x IO -27 GeV -3 

< 7 x 1CT 21 GeV" 5 


< 9 x IO" 34 GeV" 1 

< 8 x IO -27 GeV -3 

< 7 x 10- 20 GeV" 5 


< 4 x 10- 34 GeV" 1 

< 3 x 10- 27 GeV" 3 

< 2 x 10- 20 GeV" 5 


|E jm a iv(M) (fc ( ( ^ m + ifc ( %J| 
|E 3m 2 y jm (e,0) (k\% jm + ik\% jm )\ 


< io- 37 

< IO" 31 GeV" 2 

< IO" 25 GeV" 4 


< io- 38 

< IO" 32 GeV- 2 

< IO" 25 GeV" 4 


< IO" 38 

< IO" 31 GeV- 2 

< IO" 24 GeV" 4 


< io- 38 

< IO" 31 GeV" 2 

< IO- 24 GeV" 4 



TABLE III: New constraints on coefficients for Lorentz and CPT violation. The first three rows 



give constraints on general coefficients for CPT-odd operators with d = 5,7, 9. The next three 
rows display constraints within the isotropic limit. The final three rows provide the approximate 
maximal sensitivity of each source to coefficients for CPT-even operators with d = 4, 6, 8. 

opening angle 4\1> centered around q. A calculation shows that the difference ^2 — fpi in 
linear polarization at two different energies E\ and E 2 satisfies 

Sin2(^ 2 -^l) = ; sir^ec^eos^-eosfrQ 

V ' ^/(l-sin 2 2* sin 2 ^(l-sin 2 2* sin 2 <J> 2 ) 

0082(^2-^) = ; cos^+sin^cos^cos^ / g x 
Vr T ' ^(l-sin 2 2* sin 2 ^(l-sin 2 2* sin 2 * 2 ) 

This result could be used to place constraints when the observed difference tjj 2 — i>\ is small. 

An alternative and simpler strategy is again provided by considering the effective degree 
of linear polarization. In the CPT-even case, the effective degree of linear polarization 
is decreased both by polarization smearing and by the conversion from linear to elliptical 
polarization. The maximum effective degree of linear polarization is 

n efr = yjl - (l - (cos $) 2 ) sin 2 2# . (9) 

An observation of linear polarization LT then places a lower limit on U cS , which leads to 
the inequality 1 — LT 2 > (1 — (cos<3>) 2 ) sin 2 2^. A single source therefore bounds a region in 
coefficient space but cannot provide a strict constraint, as its light could be propagating in 
a normal mode with \1/ = or ir/2. In principle, combining the results of multiple sources at 
different sky locations and having different polarizations would permit a complete coverage 
of the CPT-even sector for each d, but at present the number of sources is insufficient for 
this. The above inequality can, however, be used to estimate the maximum sensitivity to 



8 



7T 

~ 2\Et 3 - W^W Y (10) 



coefficients for Lorentz violation achieved by a given source. The depletion in polarization 
is largest when the light is in an equal admixture of the two normal modes, ^ = ±7r/4. The 
effective degree of polarization then reduces to II eff = |(cos$)|, which vanishes for relative 
rotations greater than n across an energy band. We can therefore estimate the maximal 
sensitivity to the CPT-even coefficients as 

jm 

This result can be applied to the sources GRB 041219A, GRB 100826A, GRB 110301A, 
and GRB 110721 A. For Lorentz- violating operators of dimensions d — 4, 6, and 8, there 
are 10, 42, and 90 independent vacuum coefficients, respectively [10]. Table III lists the 
resulting constraints on a linear combination of these coefficients generated by each GRB. 
They represent sensitivities improved by factors of 10 to 100,000 over existing bounds on 
CPT-even violations of Lorentz invariance. 

This work was supported in part by the Department of Energy under grant DE-FG02- 
91ER40661 and by the Indiana University Center for Spacetime Symmetries. 
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